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Abstract 

We derive a total set of MHD equations in SM describing evo- 
lution of a dense plasma with neutrinos. First this is done for a hot 
pair plasma consisting from electrons and positrons, neutrinos and an- 
tineutrinos of all flavors in an isotropic medium like the early unverse 
plasma at the lepton stage. Then we find how axial vector currents 
violating parity in SM contribute to MHD for a slightly polarized 
(anisotropic) plasma where a new mechanism for the amplification 
of mean magnetic fields arises due to the collective neutrino-plasma 
interactions instead of assumed asymmetry of fluid velocity vortices 
leading to the same eff'ect of a^-dynamo. 



1 Introduction 

It is well-known that the magnetohydrodynamic (MHD) or macroscopic de- 
scription of a plasma is less detailed and much simpler than the kinetic one 
which corresponds to the microscopic description of the plasma evolution and 
therefore it is a more complicated approach. 

*Talk given by author at the International Conference on Astrophysics and High Energy 
Physics AHEP-03, Valencia, Spain, October 2003. Published in Proceedings PRHEP- 
AHEP2003/059. 
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The MHD equation system allows, in particular, to derive the Faradey 
(induction) equation for magnetic fields in the standard model of electroweak 
interactions (SM) including weak interaction terms. The main goal of this 
work is the detailed derivation of Faradey equation in SM that is important 
for the generation of primordial magnetic fields in cosmology and magnetic 
fields in a supernova protostar where powerful neutrino fluxes interact with 
the dense plasma. 

We derive the full set of MHD equations using the standard method of 
moments [1] for Relativistic Kinetic Equations (RKE) written in the coUi- 
sionless (Vlasov) approximation. There are other ways to derive MHD, e.g. 
using the Lagrangian formalism for relativistic multicomponcnt fluid [2] while 
we prefer the method of the quantum RKE for lepton plasma in SM [3, 4] 
that is more approppriate to describe both classical and spin properties of 
polarized plasmas permiated by an external magnetic field. 

Note that neutrino RKE is a useful tool to describe many phenomena in 
astrophysics and cosmology. In particular, neutrinos play the most important 
role for a supernova (SN) burst or in the lepton asymmetry formation before 
the primordial nucleosynthesis in the early universe. The usual motivation 
to use the RKE approach for neutrino propagation in a dense matter is 
stipulated by the account of neutrino collisions: within a SN neutrinosphere 
or in the hot lepton plasma of the early universe before neutrino decoupling. 

However, in addition to collision integrals there are self-consistent weak 
interaction terms in the neutrino RKE [3] that are linear over the Fermi 
constant ~ (see below section H) and analogous to the Lorentz force 
terms for charge particles in the standard Boltzman RKE which in turn are 
linear over the electric charge ^q{q = — \ e\ for electrons). 

Let us remind that these self-consistent electromagnetic fields play a very 
crucial role in the standard kinetics. In coUisionless, or Vlasov approxima- 
tion, such kinetic equations describe, e.g. thermonuclear plasmas in labora- 
tory and stars for which an energy exchange between electromagnetic waves 
propagating in plasma (=eigen modes) and charged particles (wave-particle 
interaction) proceeds faster than via the direct particle collisions (through 
particle-particle interaction) with all following issues in coUisionless plasma: 
instabilities, heating, etc. 

One expects that the presence in the neutrino RKE of the similar terms 
F™^"*^(x,t)5/('^)(k,x,t)/9k where the weak force F""^"*^ = dV""''^ / given 
by the neutrino interaction potential V^'^"'^ ~ GFne{'x.,t) is linear over ~ 
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Gp could lead for neutrinos to some analogous collective interaction effects, 
e.g. to neutrino driven streaming instability of plasma waves in an isotropic 
plasma [5] or instability of spin waves in a polarized medium [4] , and to the 
generation of magnetic fields in hot plasma of early universe [2, 6]. 

Note that in literature describing neutrino oscillation phenomena one 
neglects the changes of neutrino momentum coming from non-forward scat- 
tering, d/dk{...) —0, or refraction from density variations, d/dx.{...) =0, [7]. 
Vice versa, if one neglects neutrino masses and neutrino oscillations this 
weak force remains as a main contribution in Vlasov approximation slightly 
changing neutrino trajectories in the WKB approximation kx ^ \, where 
the momentum transfer | q |=g ~ x'^obeys the WKB condition q <^ k. 
Such small momentum transfer q = k — k' 7^ corresponds to the plasmon 
(Cerenkov) emission by a massless {K"^ = k^ — k'^ = K'^ = k^ — k'"^ = 0) neu- 
trino in medium, z/(i^) ^{K') +7*(Q), and the resulting weak force F'""^"*^ 
is the friction force acting on neutrino fluid due to this elementary process. 
Some details and the explicite form of such force (acting from plasma on 
neutrinos) in dependence on dispersion characteristics of an isotropic plasma 
can be found in the preprint [4]. 

In this work we try to find a new interesting consequence - generation of 
magnetic fields by collective neutrino interactions following from the presence 
of weak forces which are additive to the usual Lorentz force and act from 
neutrino fluid on the electron-positron plasma. For that we should derive 
the set of MHD equations and then to generalize Faradey equation in SM. 

2 Lepton MHD in Standard Model (SM) of 
electro weak interactions 

In this Section we derive MHD equations using the method [1] of moments of 
kinetic equations , or integrating RKE's over momenta, / d^p{...), J d^pp x (...) 
/ d^pCpX (...). We start from the simple case of unpolarized (isotropic) plasma 
and in the next subsection we derive MHD using RKE's in a magnetized 
plasma [4]. 
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2.1 Lepton MHD in unpolarized medium 

In an isotropic unpolarized plasma the collisionless RKE for electrons and 
positrons (e = ± | e | with upper sign for positron) derived in SM including 
weak forces takes the form [5] 

9/W(p,x,i) 9/W(p,x,i) , , , , a/W(p,x,t) 



+v X V X (j^„(x,f) -jpjx,i)) 



+ 



= 0, 



(1) 



wherej(--^'')(x,i) = K„,p„(x, i), j,„,p„(x, i)) = / d3A;(A;^/£;fe)/('^'^)(k, x, i)/(27r)3 

is the neutrino (antincutrino) four-current density; Cy^ = 2,^ ±0.5 is the vec- 
tor coupling for a = e, ij,,t neutrinos with the upper sign for the electron 
ones a = e. 

We complete the system by the neutrino (antineutrino) collisionless RKE's: 
Oji-^ + + G.V2c^^ [-V (.(^)(x, t) - .(^-)(x, t)) - 



a/(-")(k,x,t) 



dk 



0, 



9i 



^-^^ - G^V2c^^ [-V (n(^)(x,t) - n(^-)(x,t)) - 



- ^^"^(-'^^,;^"^(-'^^U nxVx(j(-)(x,t)-F)(x,t 



a/(^")(k,x,t) 



dk 



0, 
(2) 



tt)^ is 



where j(^'^")(x,t) = (ne,e-(x, t), je,e-(x, t)) =/d-V(p^/i?p)/('^'^")(p,x,t)/(2 
the electron (positron) four-current density. 

In 5-moment approximation of ideal hydrodynamics we neglect collisions 
and hence omit viscosity, heat flux terms while retaining self-consistent clec- 
troweak interactions between leptons. Thus we have to derive the particle 
density conservation (continuity) equation, the motion (Euler) equation (mo- 
mentum conservation) and energy conservation equation. 
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Continuity equations 

The weak interaction forces above have the Lorentz structure or enter 
RKE's as well as the electromagnetic Lorentz force in the third term of Eq. 
(1), F^^"'^{x,t) X (p'^/ep)9/*^")(p, X, Hence they do not contribute to 
the continuity equations [5] which take the standard form dj'^^^ / dx^ = after 
integration of RKE's (1) and (2) over momenta d^p, d^k correspondingly, 
resulting in 

dn±{x,t) ^ a[n±(x,t))V±(a;,t)] ^ ^ 
dt dx 
for charged leptons and 

dn.Mx,t) aK,..(x,t))V(--'^-)(x,t)) _ 

dt ^ dx ' ^' 

for neutrinos (antineutrinos). Here n± = n^7±, n,y^^p^ = n'^^ ^^'^^^^^^ are the 
lepton densities in the laboratory reference frame. The four-currents j^"-* = 
n'gU^^ are given by the Lorentz- invariant densities = j'^^U^"'^^ where 
U^""^^ — (7a,7aV*^"^) is the unit four-velocity of the plasma a-component, 

Uia)u(-)t^ = 1, 7a = (1 - V^)-y^ a = ±, Ua, Pa. 

Considering the particular case of the hot pair plasma = = T /j, 
where the fast e'^'-f interaction provides equilibrium leading to the zero chemi- 
cal potentials /x_ = —n+ = n = and introducing the small perturbations for 
comoving components V± = V -|- SY±, SY± -C V that means 7_ ft; 7_|_ = 7 
we get the single continuity equation for the total charged lepton density 
n' — n'_ + n'^ instead of the two ones in Eq. (3), 

d^n'{x,t) ^ a[7n-(x,t))V(x,t)] 
dt dx 

Such hydrodynamical approximation means strong correlation in a dense 
plasma between opposite charges due to which the continuous lepton medium 
becomes electroneutral conducting liquid (electrons and positrons move with 
the same velocities as a whole) resulting in electric field vanishes while mag- 
netic field exists (lepton MHD, see below Eq. (9)). Neglecting protons the 
electroneutrality condition means that the background densities n'_^Q entering 
the total ones n'_^ — n^Q -|- 5n'_^ obey the equality 

n'_o = n'j^Q = noe, 
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where in the hot plasma noe = 0.183T^. 

Note that the fluid (mean) velocity V differs from the microscopic n that 
enters the RKE of massless particles (2), | n |= 1. Of course, the Lorentz 
transformation with the unit vector Uj^"'' — {^ua^li^a^i^a) docs not change 
the value of the microscopic four-momentum — {Ek,k), El — k'^ — 0. 

Motion equations 

Multiplying the RKE (1) by the momentum p and integrating it over d^p 
with the use of the standard definitions of the fluid velocity 
V±(x, t) = n^^ / d^pvf^^\p,x,t)/(27rY, the positron ( electron) density 
^± = / '^^P/'-^^(p, X, t)/(27r)^, and the generalized momentum of the lepton 
fluid P± = w±7±V± = n^^ / d^ppf^'^\p,x,t)/{2TT)^, one obtains the Euler 
equation that coincides with Eq. (4.6) derived in [2] using another (relativis- 
tic Lagrangian) approach for multicomponent fluid, 



{dt + V± • V) P± = ^± I e I (E + [V± X B]) ^ 

n± 



V 



- V5n.„ (x, t) - ^^^1^ + V± X V X (x, t) 



(6) 



where Sui,^ = n^,^ — Up^, 6},^^ = — ha ^"^^ neutrino density and neu- 
trino 3-current density asymmetries respectively; w± = e± + p'^/n'^ is the 
Lorentz-scalar enthalpy per one particle; e±, = n'j^T± are the internal 
energy and the pressure correspondingly, T± is the Lorcntz-invariant tem- 
perature. In particular, for the Jiittner equilibrium distribution f±{p) = 
exp[(/i± — PnU^)/T±], where /x± is the Lorentz- invariant chemical potential, 
the thermodynamical characteristics are also Lorentz-invariant, 

w± = rUe^^^^P^^, e± = w±-T±, = 47rm^r|X2(me/r±) exp(//±/r±) 

J^2\Jl^el J-±) 

(7) 

For equilibrium pair plasma — T_ — T all these characteristics coincide, 

W+ = W- = We, P+ = pL = Pe, etc. 

Summing Euler equations for electrons and positrons (6) one obtains the 
motion equation (electric field and neutrino density terms do not contribute) 

d(P, + P_)^_V(p',+p'_) _ 

dt 7noe 
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-G^V2^4") [(V+ - V_) X V X 5j.„(x,i)] , (8) 

where d/dt = d/dt + V • V. 

Then we use in (8) the Maxwell equation without displacement current 
{dE/dt is omitted in MHD ), =| e | noe7(V+ - V_) = 2 | e | = 

rot B/47r where Ue = 7T^oe is the plasma density in the laboratory reference 
frame. We put also the total pressure p = p'_^ + p'_ = 2pe, the total enthalpy 
w = + = 2we introducing the total generalized momentum P = 
P+ + P_ = w-fV. 

Thus, we obtain finally the MHD motion equation for pairs generahzed 
in SM with neutrinos, 

dP Vp rotB X B GfV^ (a) r ^ ^ c-. / m 

The motion equations for neutrinos and antineutrinos are derived multi- 
plying the RKE (2) by the momentum k and integrating over (Pk, 



dt n^^ 



+ F.„, (10) 



dt rij^^ 
where the generalized momenta 

are given by the Lorentz-invariant thermodynamical functions Wv^,va — ^va,i^a'^ 
p'l, D D ; the weak forces given by 

= +G^y2ci?)[-VW^)(x,t) - M^^M + X V X 5j(^)(x,ty 

F,„ = -G^v^ci?) [-V5n(^)(x, t) - M^l^ + Vp„ X V X <5j(^)(x, t) 

(11) 

have opposite signs and, in general, depend on different fluid velocities, 
N 7^ Vp^. Here we input charged lepton density and 3-current density 
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asymmetries, dn^^^ = ri- — n^^ 5\^^^ = j*^^-* — j*^*^-* which are small in hot 
plasma. 

Since there are different fluid velocities as well as possible different ther- 
modynamical functions, w^,^ — e^^ + T,^^ — 4T^^ ^ w^^ = 4Tp^ with the 
equation of state for massless neutrinos, Pi, = e^^/S (see (7) for massless 
particles mj, = 0), we consider different motion equations for neutrinos and 
antineutrinos (10). Note that the inequality for electron neutrino species, 
Ti,^ 7^ Tp^ can arise due to beta-processes and the CC-current interaction and 
leads to a temperature difference for electron and muon (tau) neutrino com- 
ponents. For the latter (z/^, z/^ ) one expects same temperatures, however, 
we do not use this property to simplify the system (10). 

Energy equations 

Multiplying the RKE (1) by the energy Ep and integrating over (Pp one 
obtains the energy conservation law (upper sign for positrons) 



dt dx 
= ±7±n'±V±.(|e|E-Fr"'^) , (12) 

where the inner energy e±, the enthalpy w±, the pressure p± are given by 
Eq. (7); the weak force in the r.h.s. acting on charged leptons is given by 



yweak ^ GfV2J2cv 



-Vfe„.(x.t)-5«^ 



(13) 



Adding energy equations (12) and using the relation E = —V x B that 
is valid for an ideal conducting medium one gets the MHD energy equation 
for pairs generalized here in SM including weak forces, 

a[7^n^ee + 7^VV] ^ dj'y^n'weV] ^ (rot B) ■ (| e | [V x B] + F^"') 
dt dx. Att \ e \ 

(14) 

where the force F^^"*' is given by (13), n' — n'_^ + V_ = V+ a; V, 
7+ = 7_ ?5i 7 and meaning the equilibrium reached through the fast 67- 
interaction, T+ = T_ = T we put Ce = e+ = e_, = w+ = W-, P = p'+ 
in the agreement with (7). 
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Analogously multiplying the neutrino (antineutrino) RKE (2) by the en- 
ergy Ek and integrating over d^k one obtains the energy equation (upper sign 
for neutrinos) 

dt dx 

where the weak force acting on neutrinos from the pair plasma "Pvafia is given 
by Eq. (11). 

The set of MHD equations: the continuity ones (4), (5), the motion ones 
(9), (10) and the energy ones, (14), (15) is completed by the Faradey equation 
for the magnetic field B generalized in SM due to weak interactions (see next 
section, Eq. (24)). 



2.2 Lepton MHD in polarized medium 



Analogously with the case of unpolarized medium we can derive MHD equa- 
tions in the presence of a strong large-scale uniform magnetic field Bq which 
polarizes plasma populating partially the main Landau (non-degenerate) lev- 
els for free electrons and positrons. Other levels being populated by Icptons 
with opposite spin projections are degenerate (the factor Lande (^e = 2 dou- 
bles such states) and do not contribute to the medium polarization. 

The lepton density at the main Landau level in anisotropic medium is 
given by 



W) 



5n 



27r2 jQ 

(±) 



rdpft\ep), (16) 
Jo 



where in the hot plasma T ^ ji one obtains Uq ~| e | SoTln2/27r . 

Now using the electron RKE Eq. (30) from [4] we can generalize the Euler 
equation for electrons and positrons (6) for the case of a polarized medium, 

(9* + V± • V)P± = I e I (E + [V± X B]) ^ 



-VSn.A^,t) - ^^k^^ + y^ X V X 5j.„(x,i) 



^GW2 ^ ^±)^(o)di^ _ v(b(o) . 5j.„(x, t)] , (17) 
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where in a non-relativistic (NR) plasma, the relativistic polarization density 

tj Grins 

(±) _r4 I e I ^ome ± ^^(3 - t;^) 
" 3 + IStt^ io ^^^^''^ dp 

coincide with the main Landau level contributions, A^o^"* — > nQ^\ given by 
Eq. (16); B = Bq + B' is the total magnetic field. 

Adding equations (17) we obtain finally the pair motion equation in po- 
larized medium: 

dP Vp rotB X B Gf\/^ (a) r ^ c / m 

di Airrie \ e \ Anne u 



^GW2 ^ . (_) _ ^(+))b(°)div <5j.„(x, i) 



e 



-(iv(-)-iv(+))V(b(°)-<5j.,(x,t))] . 

(18) 

Note that the polarization asymmetries Uq —n[^\ "* — A'^q^'' are small 
in the hot relativistic plasma of early universe while in a degenerate electron 
gas of a magnetized supernova, T <^ /i, these asymmetries can be large since 
ni~^ =1 e I Boi^/27r^ > = (| e | BoT/27r^)e-''/^ . 

Assuming B' ^ Bq we can include the perturbative field B'(x, t) into the 
polarization terms in the last lines of Eq. (18) with the change Bq B, 
b(o) ^ b. 

In general, one can consider the limit of strong magnetic fields (or diluted 
media) for which the main Landau level is populated only. E.g. a degenerate 
electron gas obeying the condition eB > //^/2 would be fully polarized, or 
rze ~ tiq [8], that could lead to comparable contributions of pseudovector 
and vector terms in the pair motion equation (18). 

The neutrino (antineutrino) motion equations take the form which is sim- 
ilar to Eq. (10) while in a polarized medium the vector force for neutrinos 
F,^^ (11) (and similarly Fp^ for antineutrinos) is added with the additional 
axial vector force, F,,^ F^,^ + F^^^ , 

dt - + ' 

dKp^ _ _ Vp^ (,1) 
dt " n.„ +*'^<' + *'^« ■ ^'^^ 
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The latter term ( Fl^^ and similarly F^^'* with the change of common sign 
and fluid velocity Y^,^ — > Vp^), 



^(A) ^ Gfc-a 



-VMo(x, t) - ^^^pll + v., X V X 5A(x, t) 



, (20) 



depends on the spin density asymmetry 5A^(x, t), 

M„(x,() = 4->(x,«)-.4;+'(x.t) . 

J (27r)'^ £p V me me(£p + me) 



that is given by the charged lepton spin distributions S(^^(p, x, i) obeying 
the spin RKE like Eq. (12) in [4]. In NR plasma such electron spin distribu- 
tion defines the well-known hydrodynamical characteristic - magnetization 
m(x, t) =1 hb I /(o?^p/(27r)^S(")(p,x, t) k,\ \ A(x, t) which obeys the 
Bloch evolution equation and completes the system of MHD equations for 
fermions in a polarized medium [4] . 

We note here that the neutrino (antineutrino) currents j^^^p^(x, t) entering 
through weak forces the pair motion equation (18) are connected with the 
neutrino generalized momenta Kj,^ p^(x, t) via 

j..,p.(x,t) = -^K..,p.(x,t). (22) 

Note also that continuity equations (3), (4) are fulfilled in polarized medium 
[4]. We do not consider here energy equations that are easily derived from 
RKE's in polarized medium analogously to Eqs. (12-15). 



3 Faradey equation in SM 

In order to derive Faradey equation let us multiply the electron (positron) 
hydrodynamical motion equation (17) by — | e | (and + | e |) correspondingly 
and then sum them to obtain the auxiliary result for the electric field in a 
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polarized plasma: 



<T=± 



E = -^EV.xB+^^ (9,P, + V^m^V, + (V,V)P, T {V,)rN{P.l 



I ^ I i^a ^ (T=± 



(a) 



(+)^ 



dSn^^{x,t) 

at 



\/2 I e I rie — 
+ (7V(-) + 7V^))V(b-5j,„(x,t)) 

Let us stress that instead of the difference of electron and positron con- 
tributions in axial vector terms entering the pair motion equation (18) and 
given by the polarized density asymmetries ~ {tiq ^ — '^o^'') we obtained here 



the sum of them 



(-) 



+ nf,+^) that can lead to an essential effect in hot 



plasma (see below section IV) . 

Using for the last term at the first line of Eq. (23) the identity {Va)n^iPa)n- 
(Vo-V)Po- = Vo- X V X and the thermodynamics relation for the work 
dR^/dt = YAPa/dt = -p^dvjdt, (V;)„V(P.)„ = V{e^ - T^S,) + S,VT„, 
where , are the internal energy and entropy per one particle (of the 
kind a — v^, are the pressure, the volume and the tempera- 

ture correspondingly; then substituting Eq. (23) into the Maxwell equation 
9tB = —V X E we obtain the Faradey equation generalized in SM with 
neutrinos and antineutrinos: 



dfB 



VxVxB-Vx77VxB + ^ ( ^) VT^ x VS^ - 



\2e^J 



- E (^) V X {dtP, - V, X V X P,) 



^ X {dtSi.^ - V X V X + 




(24) 
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Here the equalities 6V+ + (5V_ = 0, or V+ + V_ = 2V, V+ - V_ = 25V+ = 
25V followed from the e7-equilibrium are taken into account; the magnetic 
diffusion coefficient rj — {4:7racond)''^ stems from the third term in the electric 
field (23) given by the electromagnetic collision frequency Ugrn, which enters 

the plasma conductivity acond = ^pi^^el'^e) I^T^^em with Sg, ujp = AnaUe/ Se 
being the internal energy and the plasma frequency correspondingly. In the 
non-relativistic plasma the enthalpy coincides with the internal energy, 
We ~ £e ~ "^e; while in the hot relativistic plasma — 4T, — ST. For 
the uniform conductivity the second term takes the standard form +77V^B. 

The first term in the r.h.s. (24) represents the nonlinear dynamo ef- 
fect, the third one is the Bicrmann battery effect. The fourth term can be 
neglected for small fiuctuations 5P <C P, 5V <C V. 

In an unpolarized medium we can omit all terms in last lines which are 
proportional to the axial vector coupling c^^ . The remaining standard terms 
and weak interaction vector terms (~ Cy^) reproduce the Faradey equation 
(5.7) in [2] for the ideal pair plasma (r) — 0) interacting with neutrinos 
(antineutrinos). 

The neutrino (antineutrino) currents jua,Pa entering (24) are given in Eq. 
(22) by their generalized momenta i^ua,Pa which in turn obey the motion 
equations (10), (19). 

In the next section we consider an application ^ of the generalized Faradey 
equation (24) which includes weak interaction terms violating parity to the 
actual problem of magnetic field generation in the early universe plasma. 

4 Large-scale magnetic field generation in early 
universe 

The main problem of primordial magnetic field generation that leads to a 
seed of observable galactic magnetic fields is an inconsistency of their values 
B and correlation lengths Lq obtained in the different scenarios. 

There are many ways how to generate small-scale random magnetic fields 
with large values of Brms — V< >, e.g. using some causal mechanisms 
like bubble coUisions at phase transitions, while the correlation length of such 
magnetic fields evolved (via inverse cascade) during expansion of universe 

^Results in section below were obtained together with D.D. Sokoloff in [6]. 



13 



into large-scale magnetic fields turns out to be too small at present time, 
Lq ~ tens par sees, to reach the size Lq ~ 100 hps for galactic magnetic 
field, or even more Mps) for extragalactic magnetic fields. The other 
way using inflation scenario allows, vice versa, to get large-scale (a few Mps) 
magnetic fields while their strength occurs too small for observable magnetic 
fields. 

Let us simplify the Faradey equation (24) rewriting it as a simple govern- 
ing equation for mean magnetic field evolution 

= V X aB + rjV^B , (25) 

where we omitted: dynamo term neglecting any macroscopic rotation in 
plasma of early universe, Biermann battery effect and weak interaction terms 
given by the vector coupling Cy^ because of the absence of neutrino vortic- 
ity in isotropic neutrino gas, V x j,^^ (r, t) = 0. The neutrino fiuid vorticity 
vanishes also for isotropic neutrino emission from a supernova since in the 
diffusion approximation neutrino fluxes (not a particular neutrino) propagate 
along radii, jt/„(r, i) || r even under neutrinosphere. 

In Eq. (25) we approximate the tensor coming in E from the axial 
vector force in (17) by the flrst diagonal (~ «%) term: 




(26) 



where densities Hq arc given by Eq. (16), riy/rif, = 0.5, and we assume a 
scale of neutrino fluid inhomogeneity t ~ Ag^j^j, that is small comparing with 
a large A-scale of the mean magnetic field B, Ag'^lj -C A. 

The diffusion coefficient r] ^ (47rl37 T)~^ is given by the relativistic 
plasma conductivity. 

For a small neutrino chemical potential fi^, (.uaiT) = fiu^(T)/T <^ 1, the 
neutrino asymmetry in the r.h.s. of Eq. (26) is the algebraic sum following 
the sign of the axial coupling, = ±0.5, 

^-i:41'^^||K.(r) + 6,(T)-6.(T)l. (27) 
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We stress that the Eq. (25) is the usual equation for mean magnetic 
field evolution (see e.g. [9]) with a-effect based on particle effects rather on 
the averaging of turbulent pulsations. It is well-known (see e.g. [10]) that 
Eq. (25) describes a self-excitation of a magnetic field with the spatial scale 
A ^ ri/a and the growth rate a'^/At]. 

Substituting a into A = r]/a we arrive now to the estimate 



where the neutrino mean free path luiT) = Ty^^ is given by the weak rate 



Tw = 5.54xl0-22(T/MeV)5 MeV, Ih{T) = {2H)-^ and 7i=4.46xl0-22(T/MeV)2 MeV 



is the Hubble parameter. 

If the neutrino fluid inhomogeneity scale Ag^^jj is of the order ly{TQ) ~ 
4 cm <^ Ih{To) ~ 10^ cm, we have A/lfj > 1 at the beginning of the lepton 
era (T = Tq ~ 10^ MeV). The magnetic field time evolution is given by 



where -Bmax is some seed vahic at the instant Tmax <^ Tew -^100 GeV (here 
we imbed the standard estimates of -dynamo into the context of expanding 
Universe and rely on the point-like Fermi interaction used above). 

For Afl'^^i^(T) ~ lu{T) we can estimate the index in the exponent (29) sub- 
stituting in the integrand the expansion time 

t{T) = 3.84 X 102i(T/McV)-2MeV-7\/5^ with the effective number of de- 
grees of freedom g* ~ 100 at the temperatures T > 1 GeV. Then from our es- 
timates of <y{T) and 1]{T) with the change of the variable (T/2- lO'^McV) — >■ x 
one finds the fast growth of the mean field (29) in hot plasma, a; < 1, 



given by the upper limit Xmax — 1, ^max = 20 GeV. The behaviour of ^ue{T) 
at high temperatures is unknown as well as a value of the neutrino den- 
sity asymmetry. We can state only that this value changes due to neutrino 
oscillations somewhere below T < 10 MeV not overcoming the primordial nu- 
cleosynthesis limit I C^J< 0.07 at the BBN time (T ~ 0.1 MeV, a; = 5 • 10"^) 




(28) 




(29) 




(30) 
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[11]. Nevertheless, even for \^p^ |^ 0.07 there remains an enhancement of a 
mean large-scale magnetic field -Bmax ^ T'^^J | e |<^ T'^wl I ^ I by collective 
neutrino-plasma interactions considered here, or this mechanism can be effi- 
cient and important in cosmology. This is possible for a small neutrino fluid 
inhomogeneity scale A^^^j^ <C lv{T) entering the a-parameter (26) instead of 
the assumption Ag^^j^ ~ ^i^(^) used in Eq. (30). Note that this free neutrino 
fluid parameter can vary in a wide region ^ Ag^^jj < ly{T). 

Let us remind that the inflation mechanism (with a charged scalar fleld 
fluctuations at super-horizon scales) explains the origin of mean fleld at cos- 
mological scales. However, the value of this fleld is too small for seeding 
the galactic magnetic flelds. The ampliflcation mechanism suggested in our 
paper [6] can improve this very low estimate by a substantial factor from Eq. 
(30). 

Thus, while in the temperature region Tew > T > Tq = 10^ MeV 
there are many small random magnetic fleld domains, a weak mean magnetic 
fleld turns out to be developed into the uniform global magnetic fleld. The 
global magnetic fleld can be weak enough to preserve the observed isotropy 
of cosmological model [12] while strong enough to be interesting as a seed for 
galactic magnetic flelds. This scenario was intensively discussed by experts 
in galactic magnetism [13], however until now no viable origin for the global 
magnetic fleld has been suggested. We believe that the a^-dynamo based on 
the a -effect induced by particle physics [6] solves this fundamental problem 
and opens a new and important option in galactic magnetism. 

I thank Prof. D.D. Sokoloff for the illuminating discussions and the AHEP 
group at IFIC guided by Prof. J.W.F. Valle for hospitality during part of 
this work. This work was partially supported by the program of Presidium 
RAS "Non-stationary phenomena in astronomy" . 
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